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1 Introduction 

We study in this paper the generation of a fermion number (F) asymme- 
try in the early universe due to F- violations in non-renormalizable effective 
interactions. We are motivated by the possibility of generating the cosmo- 
logical baryon number asymmetry through violation of global symmetries, 
such as global baryon or lepton numbers, by quantum gravity However, 
we will not restrict ourselves to this case. We consider the generic case of 
effective interactions due to physics at a scale A < Mp. 

We have constructed a minimal toy model that, we expect, shows the 
main features of F-generation in generic models with F-violation only in 
non-renormalizable interactions. The toy model consists of fermions b and 
d, and two scalars x cr- The i^- asymmetry is generated through out-of- 
equilibrium decays of the massive x- All other fields are taken to be massless 
for simplicity. It is assumed that other non-specified interactions maintain 
equilibrium distributions for all particles except x- 

The different scenarios of F-generation are classified by ranges of values 
of the following constants. 

• K = (Ty^/ H)t=m, the "effectiveness of decay" parameter, and 

• Kvq ~ T^^/H, the "effectiveness of the non-renormalizable interac- 
tions," where tq ~ T^^/T^. 

K is given by the ratio of the decay rate of x to the Hubble constant at 
the moment when the x bosons are becoming non-relativistic, approximately 
when the temperature is equal to their mass, M. tq is the ratio of the x decay 
width due to non-renormalizable decays through i^-violating interactions, 
r^^, to the total width, F^. For the precise definition of K and rg see 
and (p|). 

The total x decay width, F,^, is dominated by renormalizable interac- 
tions, whose coupling we caU gi. Thus F;^;(r < M) ~ glM/S-n (see ( |3.12D ). 



We call (72 the coupling constant of the non-renormalizable x interactions 
that provide F- violating x decay, namely T^^{T < M) ~ g^M^ j b^' . Taking 
g^,^ the effective number of relativistic degrees of freedom appearing in the 
Hubble expansion rate [^, -ff ~ ^^fg^T"^ /Mp, from now on to be a reasonable 
number g^, ~ 10^, we get for K and rg (see ( pSD and (|33g|)), 

K^lQ-^gl{Mp/M), r^^{g2M/giAf . (1.1) 

The other independent parameters that completely define the resulting 
F-asymmetry are, 
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* gi, the coupling of the renormahzable interactions that dominate the 
X decay rate, and 

* rj and ^, two parameters related to CP- violation in decays and anni- 
hilation processes, respectively (see below). 

We have chosen to keep these last three parameters constant at reasonable 
but arbitrary values, gi = 10~^ and rj = ^ = 5 10~^, in most of the cases 
presented below, to show the effects of the different ranges of the main 
parameters K and KrQ. In most cases the effects of changing gi, r] and 
can be easily understood. 

In order to be consistent with our Lagrangian, in which the effects of the 
physics at scale A only remain in effective interactions, we have to consider 
energies below A, thus T < A. Since also M < A, from (|]|) we deduce that 
small values oi K, K <^ 1, can only be obtained for A not much lower than 
Mp, so that M ~ A ~ Mp and K becomes not much larger than W~^gl. 
Given that tq < 1 always, in the case K <^ 1 we necessarily have Kr^ <^ 1. 
Values of K not much smaller than 1 can be obtained for any value of A. 
As we will see, the main remaining ranges for which qualitatively different 
P-asymmetry generation patterns are obtained are 

- Kro < 10-2, 

- 10-2 < Kro < 10^ 



While the usual, renormahzable scenarios of i?-violation can be classified 
with the sole parameter K ^, that mainly regulates the departure of 
the x-abundance, = n-^/n^ (n stands for number density), with respect to 
its equilibrium abundance Y^'^, it is obvious that here KrQ is also important. 
This is so because only non-renormalizable interactions generate a net F- 
abundance, Yp, whose "effectiveness" is given by Kro ~ T^^/H, as wih 
become clear below. 

2 The Model 

The schematic model we consider in this paper belongs in the so-called 
standard scenario of out-of-equilibrium decays. As is well known, in such 
scenario three elements are necessary in order to generate dynamically an 
P-asymmetry ||6|, ^]. Namely, an P-violating interaction, violation of C 
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and CP symmetries, and a departure from thermal equilibrium. The last 
ingredient is provided by the expansion of the Universe, which we assume to 
be described by the Friedmann-Robertson- Walker cosmology, in the context 
of the hot big bang model [§, P]. The other two are discussed below. 

2.1 Lagrangian 

We include two different fermion species b and d to simplify calculations by 
avoiding the presence of cumbersome interference terms that are irrelevant 
for the classification of different F-generation scenarios, which is our main 
objective. Two different bosons x with wisely chosen F-number, give 

rise to simple F-violating dimension five operators, containing xo" and two 
fermions, while allowing for F-conserving renormalizable Yukawa couplings 
of X- This is a feature we want to preserve, namely, we expect that in 
generic models the decaying boson x will have both renormalizable as well as 
non-renormalizable decays. Renormalizable Yukawa couplings of the second 
boson, (T, would only unnecessarily complicate the model, thus we do not 
include them. Therefore, trilinear couplings XX^ included in the scalar 
potential to define a non-zero F-number for a, = —2F^, while is 
defined in the Majorana type Yukawa coupling of Xi b'^dx^ , to be Fy, = 
Fd + Fb. By arbitrarily choosing F^ = Fi, = 1 one obtains F-^ = 2 and 
Ff^ = —4 from the above couplings. 

The complete Lagrangian of our toy model is, then, 

C = gi{b-dx^ + h.c.) + ^ib-dx^a + h.c.) + Vix,a) , (2.1) 

with 

V{X,(J) = -M^X^X + 93{XX(^ + X^X^<^^)+9AX^X(^^<^ ■ (2.2) 

Several comments are in order. We have assumed the existence of only one 
non-renormalizable F-violating term. It could be argued that, for example, 
quantum gravity could generate also violating renormalizable terms. Even 
though no strong argument can be given against this possibility, its effect 
would be so severe that any approximate conservation of a global number 
would be invalidated. Inclusion of other F-violating terms with dimensions 
> 5 would enormously complicate the model, and we would like to keep it as 
simple as possible. At any rate, dimension five operators will be dominant 
in most models. (Operators of dim-5 and 6 have been considered in the 
framework of SUSY GUTS in 0, |lll.) 
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Even after these considerations the term 

^{b'^dxa + h.c.) (2.3) 

should legitimately be included in (|2.1| ). However, the effects of this term 
would be the same as the 52 term in ( p. 11 ). In general 52 and g'2 will not 
be equal. If both, the (72 and (72 terms, are of the same order of magnitude, 
their combined effect would be of the same order as that of just one of them. 
If, instead, one of them dominates, the other will again not change the order 
of magnitude of their combined effect. Thus, without loss of generality we 
assume §2 <^ §2 and neglect the §2 term. 

We will consider the masses of 6, d, and a to be negligibly small with 
respect to the mass of x-, 

In spite of its simplicity, the model possesses several reaction channels, 
contributing to x-decay and inverse-decay, XX annihilation, 2^2 "point"- 
scatterings, due to contact interactions of the incoming and outgoing par- 
ticles, and 2 <-> 2 and 2 <^ 3 scatterings. Two-to-three scatterings, in 
particular, give rise to a rather large number of diagrams. 

The i^'-asymmetry is generated by decays and annihilations of x parti- 
cles, while the other processes, inverse-decays and scatterings, partially or 
completely erase any asymmetry either of dynamical origin or due to an 
asymmetric initial condition. Decays and inverse-decays are the only rel- 
evant processes at small K , K < 1. As K becomes progressively larger 
than one, other terms in the evolution equations become important. Point 
scatterings, together with inverse-decays, are the most efficient damping pro- 
cesses, while 2 <-> 3 scatterings are relevant only for extremely large values 
of K, Kro Z with our choice of values of gi, rj and ^, simply because 
only then they dominate over inverse decays before all interactions other 
than decays go out of equilibrium. 



2.2 CP violation 



The interaction lagrangian (2.1) does not give rise to a violation of CP 
symmetry. In order to model CP violations, we should introduce more 
fields and couplings into (|2.1| ). Since we want to keep our model as simple 
as possible, we shall instead resort to an explicit parametrization of the 
squared amplitudes for various processes, that respects CPT symmetry and 
unitarity but violates CP P, 0]. 



5 



The branching ratios for F-violating decays of x ^^id x particles are 
defined as, 

_ r(x ^ bda) _ _ r(x ^ Ida) 
r = , r = . (2.4) 

Notice that = T^. due to CPT invariance. When CP is violated r and f 
are different. We introduce the CP violation parameter ry, 

Wo = ]^{r -r) (2.5) 

where 

rQ = \{r + f), (2.6) 

to describe CP in decays and inverse decays. 

Besides rj, we have to consider the possibility of CP violation in scatter- 
ings. There is no CP violation in "point" -scatterings (i.e., \x'^) ~^ \bd)^ and 
crossed channels) and most of the 2 <-> 3 scatterings because of unitarity, as 
explained in appendix B. However, CP violation is possible in \hh) — > \dda) 
and its back-process, as well as in XX annihilations. Moreover, the CP- 
violations in both processes are related to each other (see appendix B) so 
that we only need to define a single parameter ^ as, 

^o-(xX bba) = ^ (cr(xx ^ bba) - a{xx ^ bba)) . (2.7) 

Therefore, unitarity relations leave only two independent CP-violation pa- 
rameters in our model, r] and ^, that we assume to be independent of x. 



3 Boltzmann Equations 

We assume that b, d, a and their antiparticles are kept in kinetic and chemical 
equilibrium with other particles through unspecified fast reactions. Thus we 
take the chemical potentials of each of these pairs of particle and antipar- 
ticle as always equal and opposite, = —^J'h■> t^d = — /^j and fia = —fi^- 
Moreover, we make the simplification of considering fif, = in view of the 
symmetry of the model under interchange of b and d. We further assume 
that only the interactions in (|2.1| ) are responsible for the evolution of these 
chemical potentials in the ranges of temperature where these interactions 
are dominant. Thus, we have. 



A(p) = fdip) = e-(^-'^'')/^, /.(p) = e-^^->^^y^ 



(3.1) 
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In the range of temperatures of interest for global-charge generation, the 
X and X particles are going out of thermal equilibrium. This provides the 
out-of-equilibrium element needed for the generation of a particle asymmetry 
[§, ^, Thus, only for large T, T » M, we can assume that fxip) will be 
close to the thermal equilibrium distribution f^{p) = e~^^/'^ . 

It is customary [e.g., fH) to scale out the effect of the expansion of the 
universe by considering ratios of the number densities rii to the number 
density of photons, Yi = nj/n^.[| Because the interactions in (2.1) insure 
that 

as it will become clear below, we will need only to solve for y^, 1^, Y^j and 
Y^. Actually, (|3.2| ) implies that the net F-number abundance Yp is just 

Yf = a (y^ - y,) ~ (y^ - y„) . (3.3) 

Since we are only interested in order-of-magnitude estimates, we will drop 
the factor of four in (^]^) from now on (which is equivalent to effectively 
taking —F^ = F^- = 1). Thus, we will write below three coupled equations 
for Yp, and y_, y+, defined as 

y± = (y^±yx)/2. (3.4) 

This choice of densities is convenient because, as we will see, the contribution 
of y_ to the evolution of y+ and Yp is small, so that we will drop all terms 
containing y_ and solve only two coupled equations. 

Let us mention some useful relations between the scaled number densities 
Yi and the chemical potentials /ij. Because, 

we have, 

Ypc^Y^-Y„ = {e-^'^/^ -e^'^/^)/2. (3.5) 

Since we are looking to produce small asymmetries, \Yp\ ~ 0(10^^'^), it is 
a good approximation to consider ^ T. Expanding the exponentials 

we get Yp ~ —j^Lfj/T and, thus, 

g±M./T ^i^Yp . (3.6) 



^We could have equally well used the entropy density s = Qtan-y instead of n^, since 
we consider the effective number of relativistic degrees of freedom (;*s to be constant over 
the range of temperatures of interest to us. 
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Similarly, we obtain e^^^/'^ = 1 =F [(Xb ~ Moreover, from ( |3.2| ) and 

the initial conditions (at high enough T) Y_ = 0, 1^ — I5 = we obtain 
Y- = Yb — Yi and, therefore, 

g±M6/T r^i^Y_ . (3.7) 



3.1 Evolution Equation for ^-Number Density 



The processes that contribute to the generation of a net number of a and to 
its damping are x decays (D) and inverse decays (ID), both renormalizable 
(R) and non-renormalizable (NR), "point" -scattering (PS) (namely scatter- 
ings due to contact interactions of the incoming and outgoing particles, such 
as ~^ bd), XX renormalizable annihilations {RA) and non-renormalizable 
ones {NRA, such as XX ~* bb'^) and their crossed channels {NRCC, such 
as XX^ ~^ bb, etc). The Boltzmann equation for is of the form [^, ^ 



dn^ 
It 



+ 3nyH 



(3.8) 



where @a are the collision terms due to the processes just listed, so that a 
stands for a = D, ID, PS, RA, NRA , NRCC. The full set of equations 
is given in appendix A. 

We will now show how to deal with the Qa terms by considering a few 
of them in detail ||2|, |3|, |^. For example, Qd+id contains the term 



e(x ^ bid2a) 



dn^dUhdUddu^ [f^ip^)\M{x 

-MPb)fd{Pd)fM\M{bda - 



X)\'^ 



(3.9) 



where M{i — > j) is the Lorentz invariant amplitude for the process \i) going 
to and dllA = {d'^PA/i'^'^)^)S{p\ — rn\)). By CPT invariance, we can 
replace \M{bda —>■ x)P by \M{x bda)\'^, and momentum conservation 
implies fb{pi)fd{P2)fa{Pa) = exp{2nb/T) exp{fi^/T)f^'^{p^), therefore 



e(x ^ 61^20-) 



dn^dUbdUddU^ 



fM\M{x^bda)\' 



(3.10) 
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where r and f are the violating branching ratios defined in section |2.2| . To 
obtain the last line in ( ^.10 ), we have used (^.61) and (^), = J f^d?p/ (27r)^, 
and the thermal average of the x-decay width T-^^, 



K2{x) 



(3.11) 



Ki{x) and K2{x) are modified Bessel functions |12|, x = M/T, and (r^)i.est 
is the x-decay width in the rest frame = 
dominated by the renormalizable interactions, 



M. We assume that is 



dUidU2 \M{x ^ 61^2) 



(3.12) 



Another term in the Boltzmann equations corresponds to the annihila- 
tion through renormalizable interactions XX ^ and its back-process, 



0(XX ^ bb) 



dU^dU^dUhdUi 



fxhmxx - 

-hfl\M{bb 



f f_ f eg j-eg 

JXJX Jx Jx 



XX)\\ 
\M{xX 



(3.13) 



Here we have used the equality \M{bb XX)P = I^(XX ^ bb)\'^, guar- 
anteed by CPT invariance. This allows us to write ( |3.13| ) in terms of the 
XX-annihilation cross section. In fact, using that /^/^ = f^f^, where 
f^ip) ~ fx'ip)^ ^'^d the definition of the thermal average of va, 



{vcr{xx 
we get 



bb)) 



dn^dU^dUhdU-, fJ^\M{xx ^ bb)\' 



0(XX ^ bb) = 



«^)2) {vaixx 



bb)) 



(3.14) 



(3.15) 



We can always relate amplitudes for 3^2 reactions to the corresponding 
amplitude for the 2^3 back-process through CPT invariance. Therefore, 
proceeding as in this example, only cross-sections for processes with two 
particles in the initial state appear in the equations. In the same way 3 — > 1 
(or 2^1) processes are related to 1 — > 3 (or 1 — > 2) decays. 
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Using the evolution equation for analogous to ( |3.8[ ) (see appendix A), 
and changing to the dimensionless variable x = M/T through the relation 
dt = {xH)-^dx, we obtain, 



dY+ 

dx 



Ox) 
xH 



[(n 



Y 



+ vroYFY^" - 2roYFY^Y^''] 



- 96ro^[2(y+-y7) + (y^-Ypy7 + y_)y_] 



-^{vcr{xx bb) 
xH 



+ va'{xx ~^ hha) + va'{xx ~^ bba)) + c.ch. 



(3.16) 



As we will see in section 4, the first line of (|3.16D contains the dominant 
terms and we will always be able to neglect the others. The 1^* and 2^^^ 
lines correspond to D and ID, the 3^^ and 4*^ to PS and the remaining 
lines to RAN, NRAN and crossed channels (c.ch.) of the NRAN. 

The prime in a' indicates that the contribution to the cross section of 
a real intermediate particle {i.e., an intermediate particle on mass-shell) 
has been removed. Scattering or annihilation processes involving a real 
intermediate particle are already taken into account by other terms in the 
Boltzmann equations, and must be subtracted to avoid double counting. For 
instance, production of a x particle near the peak of the resonance through 
inverse decay, subsequently followed by its decay. This kind of time-ordered 
sequence of processes is described by the terms in the Boltzmann equation 
corresponding to each individual process. In our example these are inverse 
decay and decay. 

The need for subtraction of pole contributions is mentioned in earlier 
papers on baryon-asymmetry generation ^ , but not described in detail in 



the literature until recently |13], after we had developed our own subtraction 
method. In order to subtract the contribution of the pole, we compute the 
Laurent expansion about = of a, or of the thermal average {av) (since 
it may be easier to compute this average with the complete cross section). 
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We then subtract the term proportional to r~^, the only negative power 
occurring in the expansion, and set = 0. We thus identify the term of 
order as the virtual intermediate particle contribution to the cross-section 
in the narrow-width approximation. This procedure is easily seen to be 
equivalent to the methods proposed in [l^ , to leading order in the coupling 
constants. As pointed out there, virtual cross-sections defined this way can 
take negative values in the region around the pole. This negative values have 
no practical effects, as we have verified numerically, when the subtracted 
cross-section appears in a damping term in the evolution equations. 

The case of XX annihilations, which are a source term and in which the 



intermediate particle is a stable fermion, is discussed in section j.t 



3.2 Evolution Equation for F-Number Density 

The evolution of n„ is due to processes that involve non-renormalizable 
interactions, the only ones that violate F-number in our model. These pro- 
cesses are non-renormalizable x decays (NRD) and inverse decays {NRID), 
"point" -scatterings (PS), xX non-renormalizable annihilations (NRAN), 
their crossed-channels (NRCC) and 2 — > 3 non-renormalizable scatterings 
of b, d and a {NRS such as bd — > bda, bd — > bda etc.). 
The evolution equation for is of the form, 

^+3^^^ = _^0^ (3.17) 

where a = NRD, NRID, PS, NRAN, NRCC, NRS. The fuh equations 
for ricr and are given in appendix A, and from them we obtain, 

^ = 2ro^[r?(y+-yf'')-ypy^«-(i + 2y^'')y_] 

_i92rox-^%^ [(y+ + y7)yi. + (2 + y+ + 3y7)y_] 

_48rox-3%l [YfY+ + (1 + 2y^'?)y„] 

Yl-{Y:^'^f-Y^]C{va'{xX^bba)) 

2(yf )2(w'(xx ^ bba)+va'{xX ^ bba)) 
+ {va'{bd bda) + va'{bd bda)) + c.ch.] } 
+more 2 <-> 3 c.ch. (3.18) 



+4 ^ 



xH 
+Yf 
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The 1^* line corresponds to the non-renormalizable decays and inverse de- 
cays, and some terms from 2 <-> 3 scatterings needed to obtain the correct 
minus sign for Y^'^, as explained below. The 2^^^ and 3'^'^ lines correspond 
to PS processes. The rest of ( 3.1^ ) corresponds to non-renormalizable an- 
nihilations and the remaining terms of 2 3 scatterings (c.ch. stands for 
crossed channels). 

If we take into account only decays and inverse decays and neglect other 
processes, the r.h.s. of ( p.l8| ) reduces to. 



(3.19) 



Notice the plus sign in front of Y^"^. It is easy to see that a minus sign is 
necessary, because no F- number should be generated if 1+ = Y^"^ , i.e., if 
there is no departure from equilibrium. The correct sign is obtained once the 



2 <-> 3 scattering terms proportional to \M'{bd — > bda)\'^ 
are included. The unitarity relation, 

J2\M{i^j)\' = J2\M{z^j)\' 



\M'{bd^bda)Y 



applied to |j) = \bda), relates scatterings and decays. 



dYi^\M{x bda)\'^ + dUidn2\M'{bd bda)^' 
dU^\M{x bda)\^ + dUidU2\M'{bd bda)\' 



(3.20) 



Thus, the above-mentioned scattering terms are proportional to the non- 
renormalizable decay rates, namely 



dUidU2 \M'{bid2 bda)\'^ - M'{bid2 bda)^ 



dli. 



|M(x* ^ Ma)p - M(x* ^ bda)Y 



(3.21) 



Therefore, the contribution of 2 ^ 3 scatterings to the source terms on the 
r.h.s. of ( |3.18D is proportional to Y^ and is given by. 



xHn. 



-{1 + 2Y.Yf) fdUidU2dU3dU^dU5 e-(-^i+^2)/r 
\M'{bid2 bsdia^)]^ - \M'{bid2 63^4^5) 
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xHn. 



xHn. 



-{1 + 2Y^Yf) [dUidU2dU3dU4dU5 e-(-^»+^4+i?5)/T 

r 

\M\h1d2 ^ 63d4(T5)|2 - \M\h1d2 ^ 63d4CT5)|'^ 

-(l + 2y_yF) [dU^dUsdUidU^ e-^*l^ 



|M(x* ^ 63(i4^5)P - |M(x* ^ 63^4(75) 



^4r/roy^<?(l + 2y_y^) 



(3.22) 



which together with ( 3.19| ) adds up to the source term on the first line of 

(HI)- 

Notice that in (3.18), besides r/ that parametrizes the violation of CP in 
decays , there appears also the parameter ^ for the CP- violation in annihi- 
lations (see equations (|2.5D and (|2] 



3.3 Parameters and Initial Conditions 

In order to solve numerically the coupled evolution equations for Yj^ and Yp 
we impose two initial conditions at a small enough value of x = xq, such 
that at xo it can be safely assumed that x and x are in equilibrium. Thus, 
we take 

y+(xo) = y7(xo), Ypix^) = Yfo , (3.23) 

where YpQ is equal to zero only for F-symmetric initial conditions. Numer- 
ically, any YpQ <C Yf(oo) ~ 10"^*^ will be equivalent to zero, and in most 
cases we take YpQ = 10"^'^ as our P-symmetric initial condition. We have 
also studied F-asymmetric initial conditions (section ^), assuming that some 
other F-generating processes have acted at earlier times. In this case, for 
large enough Kro there is an early erasure and subsequent generation of Yp 
due to the processes we consider here. 

Because our effective Lagrangian is not valid at energy scales larger than 
A, we restrict our choice of xq to values where T < A. At any rate, the 
solutions of the evolution equations are stable against variations in xq, as 
long as Xq is not too close to 1, xq < 0.1, so 

0.1 ^ xo ^ ^ ~ 30^V^ . (3.24) 

In the absence of any compelling reason to do otherwise, we take gi/ §2 ~ 1- 
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The evolution of the F-asymmetry in our model is determined by five in- 
dependent parameters (besides the initial conditions xq, Ypo). It is clear from 
the evolution equations ( |3.16| ) and (|3.1§| ) that the most suitable parameters 
to classify different scenarios for the production of Yp are: K, the "effective- 
ness of decay" parameter, Kro, the "effectiveness of non-renormalizable re- 
actions" parameter, gi the coupling constant of renormalizable interactions, 
and T] and ^, the only two independent CP-violation parameters. Therefore, 
we use 

K,ro,gi,r] and ^ , (3.25) 

as independent parameters in the evolution equations. For our numerical 
solutions we choose reasonable but arbitrary values for three of them. 



91 = 10- 



^ = 7? ~ 5 10" 



(3.26) 



and examine different ranges of values of K and KrQ. 

K is the parameter whose value determines different baryon-asymmetry 
scenarios in out of equilibrium decay models with renormalizable interac- 
tions. 



K 



(rx) 



X /rest 



1 



Hix = I) 



/g* 



Svr 



+ 192(L)3U j J ^^■'^^^ 

Here (r-^)i.est is the decay rate of non relativistic x particles, T^{T ^ M) = 
(r;^)rest (it is the decay rate in the rest frame where = M), is the 
effective number of relativistic degrees of freedom entering into the Hubble 
constant, that we take to be g^ = 100, and Mp is the Planck mass. We 
assume that the decay width of i^ dominated by the renormalizable 



decays (see ( 3.12 )), namely the first term in ( 3.27 ) is dominant, 

Mp 



K 



9l 



y/g^ M 



(3.28) 



With this assumption on F^, 
is 



the parameter tq 



T^^/T^ defined in (|2 



1 



(g2M\ 



i92.n..A; • ^'-''^ 

K and Kr^ determine the "effectiveness of the reactions" because all terms 
in {dYj^/dx) in ( |3.16| ) are proportional to K (even the last one, where 
n^{va)/H ~ Kgl/x) and all terms in (dYp/dx) in ( |3.18| ) are proportional 
to Kro (even the last one, where n-y{va)/H ~ Kr^gl/x). As we will see, 
the final value of Yf, Yp{x ^ oo), depends only on ro for K < 1 and on K 
and ro for > 10 
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Evolution of Ya. 



The evolution of y+ is mainly determined by the first term in (3.16), due to 
decays and inverse decays. We can then write, 



In first approximation, we shall set Yp = 0, y_ = in ( 3.161 ). This 
approximation will be justified below. 

The 5**^ and 6**^ lines in ( |3.16| ), correspond to XX annihilations. Non- 
renormalizable annihilation terms are always much smaller than the renor- 
malizable one, on the 5*^ line of ( |3l^ ), and can be neglected. This is due to 
the fact that {g2M / giA)'^ <C 1, and also to phase-space considerations. The 
cross-section for xX ~^ bb can be expressed as, 

{va{xx^bb)) = ^f{x) . (4.2) 

A fit to the function f{x) is shown in figure |l[ At small x <^ 1 this term is 
dominant. However, in this region closely follows Y^"^ = 1/2 — O(x^), the 
difference being of order Y-^ — Y^'' ~ O(x^). Thus, 1+ remains approximately 
constant, independently of the detailed form of the r.h.s. of ( 3.16| ). A plot of 



l+(x) for several values of K in figure |2| shows the effect of the annihilations 
term which, as argumented, is small. 

For X > 1, the annihilation rate is suppressed relative to decays, and 
the first term in ( 3.16| ) determines the evolution of 1+. It is in this region. 



in which x particles are becoming non-relativistic and Y^"^ varies rapidly, 
where large departures from equilibrium can occur. This is, then, the most 
relevant region for generating an asymmetry Yp. 

The ratio of the 5**^ line, to the 1^* term, contains (!"+ + Y^'^) < 1, 
and the factor n.y{va) / (T^) = 32giK2{x)f{x)/Ki{x)x^, where f{x), the fit 
function just mentioned, is /(x) < 10~^. Thus, for values of x not much 
smaller than 1, the 5**^ line of ( p.l6| ) is smaller than the 1^* term by a factor 
gjK2{x)/x^Ki{x) < 1. 

We also have that the ratio of the I*'* term to the one on the 3"^*^ line of 
( |3.16| ) is given by a;^/24ro 3> 1, for x > 1. A similar argument holds for the 
term on the 4**^ line. Both terms come from "point" -scattering processes. 

We have numerically checked that Y^ is always at least one order of 
magnitude smaller than Yp and 1 + , so it can be safely ignored. This can be 
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Figure 1: Cross-section for renormalizable annihilations XX 
bb function of x. 



explained by noticing that the source term in the evolution equation for YL 
(see appendix A) is proportional to roYpY^ and ro ^ 1. 

We are then left with the term proportional to rjroYp on the 2°*^ line in 
( |3.16| ), stemming from non-renormalizable decays and inverse decays. This 
term would be non-negligible only for very large values of Yp 



Yf 



1 (y+-y^^) 



y: 



eq 



(4.3) 



that would never be produced through the mechanism discussed in this 
paper. These large values of Yp could only arise from large initial values 
YpQ of the F-asymmetry, if its erasure, as discussed below, is not efficient. 
However in this case the value of Yp never significantly departs from YpQ 
and the mechanism discussed here becomes irrelevant. 

We mentioned above that we only need to find the density of x particles, 
y+, produced at x of order 1 and larger in order to estimate Yp. Let us 
consider this statement in more detail. We shall treat separately the cases 
of <C1 and K > 1. 

For K 1, all ^-number changing interactions are out of equilibrium 
at x ~ 1. Thus Y+ does not change with respect to i + (x ~ 1) ~ Yj^'^{x ~ 
1) ~ 0.5 until the x-decays occur, at a; ~ XDecay We define xoecay as the 
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Figure 2: Y^(x) for various values of K. The curves have been 
vertically diplaced for clarity. Dashed lines correspond to so- 
lutions to the full evolution equation, and solid lines to the 
simplified equation (see text). 



value of X at which x-decays enter thermal equilibrium, namely. 



and using (r^(xDecay)) ~ (rx)rost we get XDecay - (VK) ^- Since dYp/dx ~ 
y+ — Y^'^, the overabundance of x in this case is responsible for most of the 
Yp produced. 

For K > 1, instead, there is never a large overabundance of Xi because 
the reactions that change number are in equilibrium at x ~ 1. While 
X-decays are always in equilibrium, inverse decays and x-number generating 
scatterings go out of equilibrium at some point Xf at which all remaining 
X-particles decay. Therefore, because dYp/dx ^ y+ — Y^'^, the production 
of Yp happens steadily for both x 1 and x ^ 1 and, consequently, the 
production for x of 0(1) and larger gives the right order of magnitude^. 

^ To clarify this argument assume the production rate is almost constant, i.e. dYp /da; ~ 
C where C is a constant, so the increment of Yp is proportional to the increment of x, 



(r^(xDecay)) = H{xBecs.y) = H{x = 1)/xe 



Decay ' 
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An approximate analytical solution of ([4.1|), for K > 1 and all x is given 
in appendix C. We will actually only use a simpler approximation to this 
solution, valid for K > 10 and K~^/^ <C x, 



Y+{x) - Y^\x) J_ 



(4.4) 



5 Evolution of Yp: Symmetric Initial Conditions. 

In this section we give analytical and numerical estimates of the asymmetry 
Yp produced in different regimes, characterized by small or large values of 
K and KrQ. We shall be more specific about these regimes below. 



The evolution equation (3.1§|) has two source terms, corresponding to 



generation of Yp by decays and annihilations of x particles, respectively, 
and their inverse processes. We shall consider them separately, deferring 
the treatment of annihilations until the end of this section since, as we shall 
see, they turn out to give a negligible contribution compared to decays. 



5.1 Decays and Inverse Decays as a Source. 

For K ^ 1, all X- number changing processes are out of equilibrium at x ^ 1. 
1+ remains constant beyond x = 1, 

Y+{x > 1) ~ Y+{x = 1) ~ y7(x = 1) = 0.5 , 

until X s-nd X decay. This happens when the decay rate finally equals the ex- 
pansion rate of the Universe, (T-j^^Xj^ecay)) ~ -^(•^Decay); at Xj)ecay 

This overabundance of 1+ for x > 1 is responsible for most of the Yp pro- 
duced. In this case it is easy to estimate the final F-asymmetry, because 
each XX pair that decays produces a net F-number 

F^{r - r) = F^2r]ro , 

where the CP- violating parameter r] was defined in ( |2.5| ) and ( p.6p . Thus, 
we obtain 

Yp{oo) ~ F^2r]roY+{xD) = Wo , (5-1) 

AYf ~ CAx and Yf(oo) ~ Cxf, since Aa; — Xf for x/ S> 1 while the production of Yf 
that happens before a; ~ 1 can contribute at most AYf — C, which does not change the 
order of magnitude of Yf (oo). 
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since Y^{x£)) ~ 1/2, and we have effectively set = 1 for simplicity (see 
(P)). 

For K > 10, x-number changing reactions are in equilibrium at a; > 
1, therefore l+(x > 1) follows Y^'^ closely. We have to consider different 
regimes characterized by the value of KrQ. As shown below, for ~ 
and gf < 10~^ (corresponding to our choice of parameters in ( |3.26D ), the 
annihilation term is negligible as a source for Yp compared to the decay 
term for reasonable values of K and KrQ. Thus, taking into account only 
decays, we have the following results. 

* For KrQ < 10~^, we have Yf(oo) ~ rQrj. Notice that tq < 10^'^ (see 
(|3^), so this case includes in particular the case K < 1. 



* For 10 ^ < KrQ < 10 we have damping by point scatterings, 

ri[en{10'^Kro)]'^ 



Yp{oo) 



lO^K 



For 10 < KrQ < 10^ inverse decays are the dominant damping process 
and we have, 

1f(oo) 



[in{W'^Kro)]K ' 

* For very large values of KrQ, KrQ > 10^, we obtain exponential damp- 
ing of the form Yf{oo) ~ exp{— (SGvr-fTrog^)^/^}, due to 2 <-> 3 scat- 
terings. 

In the following subsections we discuss these results in more detail. 
5.2 Damping by point scatterings 

We have seen that for KrQ < 10~^ we can ignore damping terms. As 
KrQ grows larger other terms in the equation for Yp start being relevant. 
We shall take them into account in what follows, beginning with the term 
corresponding to point scattering processes. In all cases we evaluate Yp by 
quadratures, applying the saddle-point approximation where appropiate. 

We shall now assume 192KrQ > 1. The evolution equation including 
only the source term due to decays and the damping term due to point 
scatterings reads, 

dYp KUx) , , , 

— ^ = 2rjrQK-^xA{x) 
dx K2[x) 

- 192i^ro^ (y+ + y^^) Yp (5.2) 
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with, 

A(x) = y+(x) -^7(2;). 
The expression for Yp{oo) can then be written as, 

Yf{oo) = Hdu 2riKro^^u [A(u)e^] x 
Jxo K2{u) 

exp|-n-|^"dzi^^(y+(z)+y7(z))} . (5.3) 

Since we are considering decays as the source for Yp, we expect this in- 
tegral to be dominated by the contribution of the region u > 1. For this 
reason we exphcitly extracted the exponential dependence from A. Further- 
more, we have K > l/(192ro), so we can use the leading-order approximate 
expressions, 

A(u)~^-^ ; y+(u)+y_f^(M) ~2yf'?(tx) (5.4) 

Ku 

3 

=.1 ; Yl'^{u)^J^4e-'', (5.5) 



to evaluate Yp {oo). 

To leading order in K and u we then have. 



2 2 

where we defined the exponent £{u) as. 



Yp{oc) = ,l-^ J^^ .te-W , (5.6) 



£{u) = u + 96KroJ^r(^,u] . (5.7) 



Writing £ in terms of an incomplete Gamma function [12| will be useful 
below, when we consider more terms in the equation. 

The exponent (|5.7| ) is minimal at u = up given by the equation 

^/^e"^ = m^j^Kro , (5.8) 

corresponding to the epoch of "freeze-out" of the damping process [||. At 
the minimum we have, 

£{up) ~ up + l (5.9) 

£"{up) = 1 + . (5.10) 
Zup 
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We then obtain the expression, 

^/2^ r] v?p 

for the final asymmetry. 

The freeze-out epoch up is an increasing function of Kr^. We can then 
roughly approximate, 



Up ^ in [96^-Kroj (5.12) 

for large enough values of Kvq. In this way we arrive at an explicit form for 
Yf, 



Yf(oo) 



in \ 96j|Kro 



where we ignored the square root in the denominator of ( ^.11 ). 



(5.13) 



Notice that Yp(oo) in (5.13) depends on ro only through the combination 
iCro, and that the dependence is very weak. This is unlike the "free decay" 
case where we had a linear dependence on tq. The transition between the 
two regimes is, of course, not sharp. The linear dependence with ro becomes 
flatter as Kr^ grows, turning into the logarithmic form of ( |5.13| ) at about 
Kr^ ~ 10-2. 

The factor of K in the denominator represents the damping due to point 
scatterings. This suppression of the generated asymmetry as a power of K 
is similar to the damping obtained in renormalizable models 

5.3 Damping by inverse decays 

The next term we shall take into account corresponds to inverse decays. It 
is the term proportional to Yp on the first line of ( 3.18| ). The evolution 
equation now reads. 



dx K2{x) 



l92Kro^ (y+ + y^) Yf (5.14) 



We shall proceed as in the previous section. (oo) is given by (|5 
but now there is an extra term in the exponent, 

^(n)=u + ^y|(l92rQ,n)+r(^,n)) . (5.15) 
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The equation for the freeze-out point up takes the form, 

-4^e-^^ = l. (5.16) 

When li^^ < 192/ ^up, corresponding to Kvq <^ 6.5, we can neglect the 
second term in the parentheses in ( 5.16| ). In this case we recover the results 
from the previous section. 

For Kvq > 6.5 we consider the approximate equation for up (compare 
with (H)), 

V^e^^=4/f^. (5.17) 
At u = Up the exponent is minimal, taking the value, 

£{up) ~ up + 1 (5.18) 
£"iu) l-lup\ (5.19) 

where we used 5uj, > 192, and ( |5T7| ). 

Using the approximate saddle-point condition ( 5.17 ), we then obtain the 
expression for Yp{oo), 

\/27r n 1 

Yp( oo) = , for Kro > 6.5 , (5.20) 

" Kup _ 5 

2up 



to be compared with ( p. 11 



In order to obtain an explicit expression for Yp we need a solution to 
( |5.17| ). We cannot neglect the factor Up in that equation, because that would 
be inconsistent with our previous approximations. Instead, for moderate 
values of Kr^ > 6.5 we can use the matching condition with the regime of 
the previous section and write 

/ [¥ 
Up in ( 96W -^^1^0 

to obtain, 

l>(oo) = ^ , V A ■ (^-^1) 

We see that the damping effects of inverse decays are stronger than those of 
point scatterings, and that they appear at a later epoch. 
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For larger values of Kr^, Ktq ^ 6.5, we must use an iterated solution 



of (IsTtD , 



1 



2 V 2 



1 



2 V 2 



and replace it in (|5.2[l|) . 



5.4 Point scatterings revisited 

Point scattering processes give rise to another damping term in the evolution 
equation, given by the 'i"^ line of Repeating the same analysis as in 

the previous sections, we are led to the following expression for the exponent, 

f(«) = « + ^yf(i!)2r(i,„)+r(I,„)+24r(|u)) . (5,22) 

The last term in £ is new. It never dominates the exponent, however, 
and only gives small corrections to Yp{oo), not larger than 30%. We will, 
therefore, not take it into account, since we are interested only in order-of- 
magnitude estimates. 



5.5 Damping by 2^3 scatterings 

Having analyzed the effect of the damping terms corresponding to point 
scatterings and inverse decays, we are left with those related to 2 3 
scatterings. The number of diagrams for this kind of processes is large, 
making the analysis quite intrincate. 

Since the simplicity of our model does not warrant a detailed treatment 
of this problem, we shall consider only two examples which we consider 
representative of the general situation. As we shall see, this processes turn 
out to be relevant only for very large values of Ktq > 10^, which for fixed 
ro correspond to K > lO^''. 

The first reaction channel we shall consider is ba bdd (see figure |3|a) . 
Its thermally averaged cross-section, for small and large values of x, has the 
asymptotic form, 

_o 0.012 

{va) ~ <j a;4 ^ (5.23) 

G2(2.83 X 10"^ - 1.27 X m~'^£n{x^)) x < 0.1 



where we used the notation = gig'2/ ^ for brevity. 
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Adding the corresponding term to ( |5.14 ) we obtain, 



dx 



Koix) \ 



K2{x) 

l92Kro^{Y+ + Y^'^)YF 

olio 9i {v(t) , 
2 37ri^^o^ Yf{x) 



(5.24) 



For small values of x, taking into account ( ^.231 ), the last term can be written 



as, 



2"37ri^ro^ (2.83 x 10" 



9l 



1.27 X 10^^en{x^)) Yf{x) , 



(5.25) 



and turns out to be negligible compared to the term on the second line of 
( |5.24| ) for x > 10"^*^. The latter will, therefore, be the most relevant term 
in determining Yf{oo) for Ktq < 6.5, as we have seen in previous sections 
since, for reasonable values of the parameters, xq ^ 10^^'^. 

For very large values of Kvq 3> 6.5, 2 <-> 3 scatterings can be important. 
In this regime, we write 1^(00) in the form (|5.6D , with the exponent given 

by 



Six) 
£'{x) 



X + Krn / dz 



I -Km 



2zYZ''(z) + 2''3'K 



,2 Q2 



2xYr(x)+2''3n 



x2 G2 



(5.26) 
(5.27) 



The minimum of £ is found numerically to satisfy xf > 10, so we can 
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approximate the expression for £' as 



(5.28) 



Let us define x such that for x < x the first term in 8' is larger than the 
second. In this case the results of section 5.3 hold. We are interested now in 
the situation in which xp > x,so that it is 2 <-> 3 scatterings that determine 
Yf(oo). 

The value of x is given by the equation, 



1 



(5.29) 



Thus, we have, 



xigi) ~ 27.7 + In 



0.01 

11 



(51 < 0.1) 



The minimum of £ will then be at 

XF ^ {36-11 Krogj] 



(5.30) 



(5.31) 



as long as Ktq is large enough so that the consistency condition xp > x is 
satisfied. For gi =0.1, we obtain 



At X = Xp we have, 



Xp > X <^ Kro > 4 X 10^ . 



£"{xp) ^ — 

Xp 



(5.32) 



(5.33) 



The minimum of £ is therefore broad, and the steepest-descent approxima- 
tion cannot be applied. However, in view of the preceeding considerations, 
we expect Yp{oo) to be exponentially damped in this regime. 



l>(oo) 



-(367r/s:rog2)i/6 



(5.34) 



As another example of 2 ^ 3 scattering we consider the diagram in figure 
). The last term in ( |5.24 ) must now be substituted by. 



2 

vr 



■Z2 9iKrof{x)Yp{x) , 



(5.35) 
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where, 



/(^)-<! ^1032 ^ • (5.36) 



+ 0(x"*) X > 1 

As in the previous case, for small x this term can be neglected compared 
to the point-scattering term on the second line of ( ^.24 ). For x ~ 1, cross- 



section ( ^.361) is very suppressed due to the subtraction of real-intermediate- 
particle contributions. Only when x > 10 can this term be relevant. 

The analysis follows the same lines as for the previous diagram, since 
the asymptotic dependence is x~^ in both cases. This is, in fact, a general 
result; since the limit x ^ 1 corresponds to low temperatures, i.e. low initial 
energies, the cross-section in this limit is essentially determined by the final 
three-body phase-space. 

The result in this case is, 

Xi. ~ 4(i^ro5?)'/' . (5.37) 

valid for Kvq > lO"^ . Thus, we expect this 2 3 process to start being 
relevant before the previous one as Kr^ grows, and to have a somewhat 
stronger damping effect. 





Yf(oo) 


K 


ro 


XQ 


decays 


annihilations 


0.1 


10-^ 


0.1 


5.0 X IQ-^ 


3.9 X lO"^'^ 


10 


10-^ 


0.1 


4.8 X 10-* 


4.5 X 10"^" 


100 


10~^ 


0.1 


3.5 X 10"*^ 


1.9 X 10-^^ 


100 


5 X 10-^ 


0.03 


2.4 X lO-'^ 


4.9 X 10"" 


1000 


10-^ 


0.1 


1.4 X 10"^ 


8.0 X 10"^^ 


1000 


5 X 10"^ 


0.03 


2.0 X lO"'^ 


1.9 X 10"" 



Table 1: Numerically obtained values of Yf(oo) for both, x decays/inverse 
decays and XX annihilations as source terms. Fixed parameters are specified 
in { ^1^ . 
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5.6 Annihilations and their back processes 



The analysis of the term on the 4*'^ hne of ( 3.18 ), due to XX annihilations 
and their back processes, as a source of F-number generation can be carried 
out along lines similar to those of the previous sections. We shall only quote 
numerical results here. As expected on physical grounds, and can be seen 
from table |l], this term is quantitatively much less important than decays 
and inverse decays. 




Figure 4: (a) xX annihilation diagram, (b) x decay followed by point scat- 
tering. 



A Feynman diagram for xX annihilation is shown in figure 1(a). The 
intermediate fermion can be on mass-shell, which makes the tree-level cross- 
section singular. Taking as a variable (see fig. @(a)), at fixed center-of- 
mass energy ^/s, the kinematical region where the intermediate particle can 
be on-shell is given by 

^-^^^ < pI < ^ + . (5.38) 



Notice that this interval is completely contained within the kinematical do- 
main < P4 < -v/s/2. 

The singularity of the amplitude at tree level corresponds to a space- 
time-ordered sequence of x decay and point scattering, illustrated in figure 
^(b), in which the intermediate particle propagates over macroscopic dis- 
tances. The fact that singularities in the physical region represent ordered 



sequences of processes has been proved in general in [14|. Clearly, the decay 
width of x has to be included in the propagator in order to obtain a finite 
cross-section. 

The differential cross-section da/dp^ for xX annihilation is shown in fig- 
ure H, for = 4 X 10~'^M which corresponds to the values of parameters 
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0.001 r 

0.00001 r 



-L- -LU 0.2 0.4 0.6 0.8 1 

Figure 5: XX annihilation differential cross section for s = 5, = 0.0004 
and M = 1 (G = (71(72/ A). The vertical lines show the boundaries of the 
kinematical domain of real-intermediate-particle exchange. 




327i^sr da 




■ dpi 










Pi 



0.305 0.31 0.315 0.32 



Figure 6: Solid line: same as previous figure. Dashed line: contribution of 
real particle exchange only, computed at F^ = 0. Dot-dashed line: contri- 
bution of virtual particle exchange only, at F^ = 0. 
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given in ( |3.26| ). It is also shown in figure ^ together with its value at zero 
width, and the differential cross-section for exchange of a real intermediate 
particle, also computed at Ty. = 0. The former is singular at the boundaries 
of interval ( 5.38D , and the latter is zero outside that interval. 

In order to estimate the F-number generated by this process, we set 
the subtracted differential cross-section da'/dp^ {i.e., the virtual-particle 
exchange differential cross-section) to be equal to dcj/dp\, with the above- 
mentioned value for T^, outside the kinematical limits for real-intermediate- 
particle exchange and zero inside that interval. 

The evolution equation for Yp, with XX annihilations as the only source 
term, can then be numerically solved with the results shown in table |l[ This 
gives an estimate of the importance of this term relative to the other source 
term. The i^-number generated in this case is negligible compared to that 
arising from decays and inverse decays, except for very large values of K and 
very small values of rg. In this last case both the decays- and annihilations- 
generated F-numbers are small themselves, since the system never departs 
much from equilibrium. 



6 Asymmetric initial conditions. 

We consider now the case in which the initial value for the density Yp is 
different from zero. To be concrete, we take Ypo > 0. It is clear from the 
foregoing analysis that the processes which will be relevant to the erasure 
of YpQ are point scatterings. Therefore, we have, 

(oo) = Ypo exp |-192i^ro J^dz 1 (y+{z) + (z)) | + Yp''^{oo), 

(6.1) 

where Yf'^{oo) refers to the value of Yp{oo) obtained in the previous sec- 
tions with Ypo = 0. 

We can easily evaluate the exponent as. 



nzl(y+(z)+y7(z)):^- (6.2) 

Jxo 2 ^ ' Xo 



to write, 



y^(oo) y^oexp ( -i^^) +Y^^^{oc). (6.3) 



XQ 

As mentioned in section |3.3| , the value of xq cannot be chosen arbitrarily 
small, since lagrangian (2.1) is not applicable at energies higher than 0{A). 
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More precisely, tree-level diagrams such as non-renormalizable decays and 
point scatterings violate the unitarity bound at an energy scale of order 
~ A/(^2- The appearance of 52 here should not be surprising, since A enters 
£ only in the combination 52/ A. 




Figure 7: Evolution of Yp from a non-zero initial value (solid 
lines). The dashed line shows Yf{x) for K = 5 in the symmetric 
case Ypo = 0. xq = 0.04. 



Thus, the minimal value of xq we can choose is, 



xq ~ ~ 30fi(i^, (6.4) 



where we used ( 3.24 ) in the last approximate equality. With this value of 



xo we have, finally. 



Yf{oc) ^ y^oexp ( "^^^ ^) +y7'"(oo). (6.5) 



gi 

Depending on the values of the parameters K, tq, gi and r] (since Y^'^™(oo) oc 
r]), and of YpQ, one of the two terms will dominate. Only when the first one 
is much smaller than the second will this model display dynamical, initial- 
condition-independent generation of global charge Yp. This is illustrated in 
figure |7[ for if = 30 and 5. In the first case the initial condition is completely 
erased, and only partially in the second. 
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7 Final Remarks 

We considered a schematic model for dynamical generation of a global 
charge. The interactions giving rise to the violation of F-number are given 
by an operator of dimension 5, whereas the renormalizable piece of the 
Lagrangian is F-symmetric. Our model is a special case of the standard 
out-of-equilibrium decays scenario (2|, P, with the particularity that the 
heavy x boson can decay through two channels only one of which, the one 
mediated by the effective interactions, violates conservation of F-number. 

We introduced CP violation through a parametrization of the matrix el- 
ements for X decays and XX annihilations consistent with unitarity and CPT 
(to the order considered here in the coupling constants) |3|, Q . We have, then, 
two independent CP violation parameters, rj and ^, which represent the net 
F-number generated in decay and annihilation reactions, respectively. All 
other possible violations of CP are related to r] and ^ through unitarity and 
CPT. We made one additional assumption, that ^ is constant. Genuine 
CP violation would require explicit modelling ||2|, ^, with the introduction 
of more fields and coupling constants. We preferred not to do so, to keep 
our model as simple as possible. 

The model assumes the standard Friedmann-Robertson- Walker cosmol- 
ogy P, ^, ^, Q , which enters the evolution equations explicitly through the 
Hubble parameter H. Furthermore, other species and interactions are as- 
sumed to exist, which do not change F-number, give a, ^ 0(100), and 
maintain light degrees of freedom in local thermal equilibrium. In order to 
establish the evolution equations and make them tractable we made two 
approximations that are known to be valid in more realistic models such as 
GUT-based baryogenesis models |^, ^. These are the neglect of degen- 
eracy factors and the use of Maxwell-Boltzmann equilibrium distributions. 
We also neglected processes of order higher than the first in gi and g2- In 
view of the negligibly small effects of 2 ^ 3 scatterings for all reasonable 
values of parameters, this approximation is justified. 

The model has five parameters, which can be chosen as explained in 



section 3.3. We kept gi, rj and ^ fixed (see 3.26), and studied the final value 
of the global charge for broad ranges of values of K and tq. Since we set 
the coupling constant for non-renormalizable interactions g2 = 0.1 = gi, the 
condition M < A leads to an upper bound for tq, tq ^ 5 x 10~^. 

Besides those five parameters, initial conditions xq, YpQ must be chosen. 
In the symmetric case YpQ <C 1^(00) the final value of Yf(oo) is insensitive 
to the previous value of xq. We have used xq ~ M/A in this case. For 
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asymmetric initial conditions the choice of xq is more dehcate. Whether 
large initial asymmetries are completely erased or not depends on the value 
of Ypo, K and ro, but also on xq. Due to the nature of our model, we 
cannot approximate xq — as in other simple models since some of the 
terms in the evolution equation (notably point scatterings) are singular at 
x = 0. Thus, we had to determine the minimal value xq can take. This is the 
value at which non-renormalizable processes such as point-scatterings start 
violating unitarity bounds, xq ~ g2M/A.. For x much less than this value, 
the description of F violating interactions as effective operators ceases to be 
valid. 

The case of symmetric initial conditions is the more interesting, since 
it displays purely dynamical generation of global charge. The appearance 
of an additional parameter vq in our model, besides K, which determines 
the effectiveness of damping, introduces some differences with respect to 
renormalizable models. As shown in section ^, damping is effective only 
when KrQ > 10~^, so even for relatively large values of K we can have 
undamped generation of F-charge. The largest Yf{oo) attainable in this 
case is r/ro. 

Also as a result of having a branching ratio ro, exponential damping due 
to 2 — > 3 scatterings is effectively postponed to very large values of K, since 
KrQ itself must be large for these processes to be in thermal equilibrium. 
We mentioned in section ^ that the number of 2 — > 3 scattering channels 
is rather large in our model, and considered two specific examples which 
showed that these processes are irrelevant for K <, 10^. (Notice, however, 
that they are important in deriving the Boltzmann equations, as explained 
in section § |^] ) . 

An important class of processes not present in renormalizable models is 
two-body point scatterings. (These have been considered in the framework 
of SUSY GUTs in ||l^, |l^.) These scatterings turn out to be the most ef- 
ficient damping process at low x. Thus, they are crucial for erasing initial 
asymmetries. For symmetric initial conditions, 2^2 scatterings give alge- 
braic damping (~ ^/K) for KrQ > 10~^. Inverse decays are also important 
for KrQ ^ 7. They are in equilibrium for x ~ 1, for large enough KrQ, and 
also produce 1/K damping of the generated Yp. 

Besides x decays and inverse decays, XX annihilations and their back 
processes are also a source of i^-generation. Since CP violation in this case 
is independent of that in decays, it is always possible to choose parameters 
in such a way that annihilations be the dominant source term (by setting 
^ ^ r/). In the absence of explicit modelling of CP violation this situation 
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cannot be ruled out. We considered this possibility artificial, and used ^ ~ r/ 
in our computations in order to be able to compare the relative importance 



of both source terms. As shown in section 5.6, the annihilations term is 
negligible compared to decays over a wide range of values of K and tq. 
When K is large and tq small enough to block damping {KrQ < 10~^), the 
i^-number generated by annihilations is comparable to that generated by 
decays, typically one order of magnitude smaller. 
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A Evolution equations 

Here, we give the full evolution equations of n^,nfj.,ni, and show that the 
i^'-number asymmetry only depends on — Y^. We define 



where aia2 ■ ■ •, • • • label different particles, and Q' is given by a similar 
expression where M' replaces M (namely the amplitude with on-mass shell 
propagators substracted) . The full evolution equations is, 

+ Huy. = -@d+id{x ^ bd) - @NRD+NRiDix bda) 

-Qpsi^X ^ bd) - 2Qps(bx ^ da) - 2@an{xX ^ bb) 
-"^^NRAN+NRCcixX ^ blu) - 2e'jvijAiV+JViicc(xX ^ bla) 



where the factor of 2 in PS and AN [R or NR) terms accounts for similar 
processes with b and d exchanged. The evolution equation is obtained 
by replacing in (A. 2) all particles by their antiparticles. 
The full evolution equation is, 

^ + Hnb = +OD+iDix^bd) + QNRD+NRiDix^bda) + QANixX^bb) 



+Ops{dx ^ ba) - Qpsibx ^ da) + @ps{^X ^ bd) 

+&NRAN+NRCC{XX ^ bba) + &NRAN+NRCC{XX ^ bla) 

-@NRCcibb ^ axx) - &NRCc{bb ^ axx) 

+&NRCc{(^X ^ bbx) + @NRCc{(^X ^ bbx) 

+&NRCc{^X ^ bbx) + &NRCc{^X ^ bbx) 

—QNRs{bb ^ dda) — @NRs{bb dda) + @NRs{dd -fr^ bba) 

+QNRs{dd bba) + @NRsida ^ bbd) + QNRsida bbd) 

+QNRs{da ^ bbd) + @NRsida ^ bbd). (A.3) 




X [/ai/aa " " " \M{aia2 > f^lp2 ' ' 

-fpJp^...\M{PiP2---^aio^2---)\\ 
= -e(/?i/32 • • • ^ aia2 • • •)> (■ 



(A.l) 



~2@NRCc{crXX ^ bb) - 2G nrcc{^XX ^ bb) 
-2@NRCcibxX ^ ba) - 2@MRcc{bxX ^ ba) 
-2Qnrcc^XX ^ ba) - 2QNRCcibxX ^ ba), 



(A.2) 
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Similar expressions hold for (switching b and d) , and for , nj (replacing 
the particles by their antiparticles. Using (A. 2) and (A. 3) one can show 
rib — iT'l ^'^d n^ — n^, evolve equally, i.e. 



d{nb - ni) 
dt 



+ H(nb- ni) 



d{nd - nj) 
dt 



+ H{nd- nj) 



+&d+id{x ^ bd) - Qd+id{x <^ bd) 
+&nrd+nrid{x ^ bda) - Qnrd+nrid{,X ^ dd^) 
+@ps{^X ^ bd) - Qpsi^X ^ bd) 
+2Qps(bx ^ da) - 2eps{bx ^ da) 



j^{n^-n^) + H{n^-n^] 



(A.4) 



Hence 



d{Yb-Y^) 



d{Y,-Y^) 



d{Yd-Ya) 



(A 5) 

dt dt dt ' ^ ' ' 

and since Yp = 4{Y^ - Y^) + 2(Y^ - Y^) + (Y^ - Yj,) + {Y^ - Y^), Yp is given 
by Y^ - Y^. 

The full evolution equation is 



dria 
~df 



+ Hua = +'dNRD+NRiD{x ^ bda) + 2@ps{bx ^ da) + @ps{bd ^ ax) 

+'^®'nran+nrcc(.XX ^ bba) + 2@NRCc(.bx ^ bax) 
+2QNRCc{bx ^ bax) + '^^NRCc{bx ^ bax) 
+2QNRCc{bx ^ bax) - &NRCcicrX ^ bbx) 
-2Qnrcc{(^X ^ bbx) + 2Q NRCc{bb ^ axx) 
-2&NRCc{ba ^ bxx) - 2eNRCc{ba ^ bxx.) 
+Q'NRsibd dba) + &'NRsibd bda) 
+2eNRs{bd^ bda) + 2eNRs{bb ^ dda) 
-'i^NRsibcr ^ bdd) - 2eNRsiba ^ bdd). (A.6) 

Changing all particles by their antiparticles in (A.6) one obtains the equation 
for n^. 

The evolution equation of Y- is, 



dY. 
dx 



xH 
-24tqx 



[Y_{1 + 2r^«) + roYpYl\l - 2roY_)] 

[y-(i + 2yf«) + yj.y+)] 
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_96rox-^%^ [y_(2 + Y+ + 2^!'') + YFiY+ + yf«))] .(A.7) 
Id 

B CP violation 

All the processes we have considered in this paper arc upto the order ((72/A)^. 
In this Appendix, we will use unitarity to check the possibility of CP vio- 
lating differences [|M(a 6)^ — \M{b a)p] upto the order (52/A)^. 

In order to have CP violation in a process \a) \b), there must be other 
final states besides \b) for |a) to go to and other initial state besides \a) that 
can go to 1 6). Namely, unitarity imposes the following relations, 

^\M{a^j)f = Y,ma^j)f (B.l) 
j j 

j2\M{i^b)\' = j2mi^-b)\' (B.2) 

i i 

where a,b,i,j denote the CP conjugates of a,b,i,j respectively. Therefore, 
if 1 6) is the unique state of \a) or \a) the unique initial state for |6), then 
either (B.l) or (B.2) imply 

|M(a^6)P = \M{a^b)f , 

which is equivalent to the requirement of CP conservation in the process. 
Using this condition, we find most processes do not violate CP upto the 
order of ((72/A)^. 

Let us look at some examples: 

1. Consider \bx) as initial state. By (B.l), we get 



/ 



dUddU^\M{bx da)f + dU^dU^dU^\M{bx ^ bax)f 

+dU-^dU^dU^\M{bx ^ bax)f 

dU^dU^\M{bx da)\'^ + dnhdn^dn^\M{bx ^ bax)\'^ 

+dUbdU^dUy^\M{bx bax)f] (B.3) 

that may seem to allow CP violation. But with \bax) and \bax) as final 
states by (B.2) we get 

\M{bx^bax)f = \M{bx^bax)f, (B.4) 
\M{bx^bax)\'' = \M{bx^bax)\'' (B.5) 
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and using (B.3), we have 



\M{bx^da)f = \Mibx^da)\' 



(B.6) 



Therefore there is no CP violation in any of the three processes in (B.3). 

2. With the assumption of no asymmetry of b and d in all processes upto 
the order of i.e. 



\Miax ^ bbx)f = max ^ ddx)\' 



(B.7) 



we find again no CP violation in processes that involve the state \crX). In 
this case, the relation (B.l) is 



/ 



/ 



dUbdUdlMiax ^ M)P + dUbdnidU^\M{ax ^ bbx)\'^ 

+dnddnsdn^\M{ax ^ ddx)f 

dUidn^\M{ax ^ bd)\^ + dUi,dUidU^\M{ax 
+dUddUjdU^\M{ax ^ ddx)\^ 
Now, consider \bbx) as the final state then (B.2) gives 

\M{ax ^ bbx)\^ = \M{ax ^ bbx)\^ = \M{ax ^ bbx)\^ 
and with \ax) as initial state, (B.l) and (B.7) yield 
\M{ax ^ bbx)\^ = \M{ax ^ ddx)f = \M{ax ^ bbx)^ = \M{ax 
Then, combining (B.9) and (B.IO), one obtain 

\M{ax^bbx)f = \Miax^bbx)f, 
and (B.ll) together with (B.8) give 

\M{ax^bd)\^ = \M{ax^bd)\^ 



(B.8) 



(B.9) 



ddxt 
(B.IO) 



(B.ll) 



(B.12) 



To sum up, we find that only the processes |x) — \bd), \x) — \bda), 
\bd) ^ \bda), \bb) ^ \dda), \dd) ^ \xx) ^ \dda), \xx) ^ IM, 

and their CP conjugates violate CP. We have related the CP violation in 
Ix) \bd) and \x) — > \hda) and their CP conjugates by using the branching 
ratios r and f. The CP violation in |x) — \bda) and the violation in 
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\bd) — > \bda) are related in ( 3.20| ). With \ba) as the final state in (B.2), we 
obtain 

dU^dU^\M'{xX bba)\'^ + dUddU^\M{dd bba)]"^ 
dn^dU^\M'{xx bba)\'^ + dUddU^\M{dd bba)]"^] ,(B.13) 
namely, 

"dU^dU^ [\M'{xX ^ bba)\^ - \M'{xX ^ bba)\'' 
dUddU^ \\M{dd - \M{dd bba)\'^ 

= 2CjdU^dU^\M{xx^bba)\^ (B.14) 

where ^ is defined in (|2.7D . We have a similar equation for \dda) as final 
state. Note that in ( pl^ ) we take \M{xx ^ dda)\^ = \M{xx ^ bba)\^, 
since we assume no asymmetry of b and d in all the processes upto order 
(92 /Ay. 



C Approximate solution for 

The equation for A(x) = 1+ — Y/^'^, neglecting terms proportional to Yb and 
y_, is given by, 

dA dV"^ 

^ = -Kg{x)A(x) - ^ , (C.l) 
dx dx 

where we used the notation g{x) = xKi{x) / K2{x). Defining ^k{x) = 
KA{x/K), we get the equation. 



dAxix) 
dx 



X 



dY^ 
dx \K 



(C.2) 



with the initial condition Ak{Kxo) = KA{xo) = 0. 

For very large K, we can expand the coefficient functions g{x/K) and 
dY^'^ {x / K) / dx in this equations about x/K = (i.e. for K ^ x), solve the 
resulting equation for Ax, and transform back to A to obtain. 



w X , / o ^ ,2 5 Kx^, n ^ ,2 5 Kx^„,\ 

AW = l/8^.SF,(-.j,— )-4.fl(3,3.^)jx 

exp(-Kx^/6) 



(C.3) 
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where iFi is a confluent hyper geometric function [12|. 
For Kx^ > 1 (i.e. K~^/^ < a; < JC) we obtain, 

" 2K-X 

which coincides with the result (see ||], eqn. (6.29)), 

A(x) = y7(x)/(i^x) 



(C.4) 



(C.5) 



if we approximate, Y^'^ ~ 1/2, valid for x < 1. Notice that in our model x 
is a scalar, so Y^'^{x ^ 1) = 1/2, instead of 1. This suggests that we take, 
for large K, 



1, 



A{x) = -y7(x) (^xSFi( 
exp{-Kx^/6). 



2 5 Kx^ 



3' 3' 6 



i-t^Uo' 



,2 5 



3' 3' 6 



(C.6) 



This expression is a very accurate representation of A for large values 
of K, which already at X = 10 departs by about 5% from the numerical 
solution of equation ( C.l ) at x ~ 1, the region where the error is largest. 



